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BTATUS  OF  THE  NAVAL  OCEAN  SYSTEMS  CENTER' B 
LONG  WAVE  PROPAGATION  CAPABILITY 


Jerry  A.  Ferguson 
Naval  Ocean  System*  Center 
San  Diego,  CA  92152-5000 


INTRODUCTION 

This  paper  describee  a new  software  version  of  tha  Naval  Ocean  Systems 
Center's  Long  Wavelength  Propagation  Capability  (LWPC) . This  program  applies  the 
concept  of  the  earth- ionoephe re  waveguide  to  vary  low  frequency  (vlf:  10  kHz  to  30 
kHz)  and  low  frequency  (1ft  30  kHz  to  60  kHz)  radio  propagation.  The  propagation 
capability  has  been  incorporated  into  a FORTRAN  program  whiah  seta  up  calculation  of 
mode  parameters  along  arbitrary  propagation  paths  or  for  user  defined  operating 
areas.  In  the  latter  aase,  the  paths  are  automatically  selected  to  account  for  areas 
of  low  aonduativity.  Previous  dependence  on  VMS  has  buen  eliminated.  Substantial 
improvements  have  been  made  in  storage  and  identification  of  the  program  outputs.  A 
set  of  programs  have  been  developed  for  graphically  displaying  the  output  from  the 
program. 

Tha  first  version  of  the  LWPC  was  a collection  of  separate  FORTRAN  programs 
linked  together  in  operation  by  a command  procedure  written  in  Digital  Equipment 
Corporation's  VMS  operating  syetem  oommand  language  (Ferguson  and  Snyder,  19B9n,b). 
The  initial  version,  designated  LWPC-0  in  this  papsr,  was  originally  a tsst  of  a 
concept,  making  full  integration  of  the  constituent  programs  premature.  The  core  of 
the  new  LWPC,  designated  LWPC-1,  is  embodied  in  a single  FORTRAN  program  named  LWPM 
(short  for  Long  Wavelength  Propagation  Model) . This  program  implements  all  of  the 
features  of  the  original  VMS  version.  In  addition,  LWPC-1  Includes  a number  of 
auxiliary  FORTRAN  programs  which  provide  summaries  of  the  contents*  of  files  and  can 
be  used  to  update  and  graphically  display  the  contents  of  its  output  files. 

The  propagation  model  implemented  in  the  LWPC-1,  which  is  identical  to  that  in 
LWPC-0,  treats  the  space  between  the  earth's  surface  and  the  lower  ionosphere  as  a 
waveguide.  The  upper  boundary  of  this  waveguide  is  characterized  by  an  losotropic 
conductivity  which  increases  exponentially  with  height.  The  exponential  increase  is 
defined  by  a log-linear  slope  and  a reference  height.  The  LWPC  incorporates  a model 
which  defines  an  average  value  of  the  slope  and  reference  height  which  depends  on 
frequency  and  diurnal  condition.  The  height  of  the  nighttime  ionosphere  over  the 
polar  caps  is  lowsr  than  it  is  at  middls  and  squatorial  latitudes.  This  modsl  has 
bsen  compared  extensively  with  available  measurements.  The  lowsr  boundary  of  the 
LWPC  is  based  on  the  Westinghouse  Geophysics  Laboratory  conductivity  map  (Morgan, 
1965) . Both  boundaries  of  the  earth-ionosphere  waveguide  are  known  to  have  possibly 
large  uncertainties  in  aertain  regions  and  seasons. 

A typical  usage  of  the  LWPC  is  to  generate  geographical  maps  of  signal  avail- 
ability for  coverage  analysis.  The  program  makes  it  easy  to  set  up  such  problems  by 
automating  most  of  the  required  steps.  The  user  specifies  the  transmitter  location 
and  frequency,  tha  orientation  of  the  transmitting  and  receiver  antenna  and  the 
boundaries  of  the  operating  area.  The  program  automatically  selects  paths  using  a 
coarse  resolution  of  15  degress  to  ansurs  that  the  operating  area  is  fully  covered 
and  then  using  a tina  resolution  of  3 degrees  to  ensure  that  all  significant  law 
conductivity  areas  are  included.  The  diurnal  conditions  and  other  relevant  geophysi- 
cal parameters  are  then  determined  along  each  path.  After  the  mode  parameters  along 
sash  path  are  determined,  the  signal  strength  along  each  path  is  computed  and  is 
then  mapped  onto  a grid  overlaying  the  operating  area.  The  data  in  the  grid  include 
a user  specified  (constant)  standard  deviation  of  the  signal  throughout  the  grid. 
Tha  signal  strengths  in  this  grid  can  be  displayed  by  other  programs  to  be  described 
later. 
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LWPC-i  incorporates  a number  of  improvement*  over  LWPC-Q  including  reduced 
dependence  on  the  operating  system,  better  integration  of  the  intermediate  data 
storage  and  modification  of  the  input  to  use  control  strings  coupled  with  input 
parameters.  An  option  for  previewing  the  propagation  paths  to  be  processed  by  the 
program  is  included.  A map  with  major  land  masses  graphically  filled  is  included  in 
this  preview  option  as  well  as  in  tha  geographical  coverage  displays.  The  main 
program,  LWPM,  incorporates  asssntial  elements  of  several  of  tha  programs  used  in 
LWPC-Oi  PRES EG  (which  sets  up  the  paths),  MODEFHDR  (which  determines  the  initial 
solutions  for  each  aat  of  segments) , tha  Segmented  Waveguide  program  (which  deter- 
mines the  solutions  for  segments  of  similar  type) , FASTMC  (which  doss  tha  mods 
summation  along  tha  path)  end  OPARSA  (which  sets  up  the  data  for  plotting  coverage 
araa  maps) . While  LWPC-0  generates  a log  fils  and  a mode  parameter  data  file  for 
each  path  which  is  processed,  lwpc-1  generates  only  one  log  fils  and  ona  data  file, 
making  fila  management  easier. 

Three  noiee  models  are  currently  available!  'ITSN'  names  the  ITS  noiee  model 
of  Zeohariaen  and  Jonee  (1970) , 'NTIA'  names  tha  noise  model  of  Spaulding  and  Wash- 
burn ( 19B5)  and  'DECO*  names  the  DECO  Waetinghousa  thunderstorm-baaed  noiee  model 
(Maxwell  «t  al.,  1970).  Tha  NTIA  model  is  identical  to  the  new  CCIR  noiee  model 
(CCIR,  1986).  The  model  named  ITSN  is  the  implementation  of  CCIR  322  (CCIR,  1963) 
which  mapa  the  basic  noise  map  parameter.  Fam,  in  Universal  Time.  The  model  named 
K”IA  ie  the  new  noise  modal  developed  using  additional  measurements  and  has  einoa 
become  tha  new  CCIR  model  described  in  CCIR  Report  322-3  (CCIR,  1986) . These  two 
models  are  baaed  on  surface  mappings  of  measurements  at  a limited  number  of  eitas. 
The  DECO  modal  of  atmoepharia  noise  ia  limited  to  the  VLF  regime  (10  to  30  kHs)  and 
ueaa  a data  baas  of  thunderstorms  and  doaa  propagation  calculations  from  theaa 
thunderstorms  to  tha  receiver  sites.  This  modal  was  calibrated  using  the  earns  meas- 
urement aitaa  which  were  used  to  develop  the  other  two  models.  Because  the  modal 
ueea  propagation  calculations  instead  of  polynomials,  it  ia  a much  longar  running 
modal  and  is  used  infrequently. 


PROGRAM  CONTROL 

The  LWPC-1  uses  character  atringa  for  program  control  and  to  specify  options. 
The  control  strings  have  the  same  meaning  and  use  amongst  all  of  the  programs.  On 
input,  moat  control  strings  may  be  abreviated.  For  example,  the  string  "TX-DATA"  can 
be  entered  in  upper  or  lower  case  and  can  be  shortened  to  "TX-D".  If  a control 
string  is  shortened  too  much,  it  will  not  be  recognized  and  execution  will  etop.  To 
make  it  easier  to  read  control  atringa  which  are  composed  of  more  than  one  word, 
dashes  are  used  to  separate  the  words,  such  as  the  abova  mantioned  "TX-DATA".  A 
blank  in  the  first  column  of  a line  of  data  causes  the  string  to  be  treated  as  a 
comment  line,  allowing  the  user  to  annotate  run  streams  for  documentation  and  to 
provide  prompts  for  editing. 

After  the  necessary  control  strings  and  the.ir  associated  data  are  specified,  a 
specific  control  string  named  "START"  is  used  to  initiate  the  calculations.  The 
program  LWPM  first  creates  a statue  file.  This  file  is  named  with  the  extension  STA 
and  contains  a lint  of  the  bearing  angles  and  lengths  of  all  of  the  paths.  As  the 
calculation  of  mode  parameters  along  aach  path  is  completed,  the  parameters  are 
written  to  a file  named  with  the  extension  MDS  and  the  corresponding  entry  in  the 
statue  file  ie  updated  with  the  date  and  time  of  completion  and  the  CPU  time  used. 
Calculations  for  successive  paths  continues  automatically.  If  the  computer  run 
aborts  for  soma  reason,  then  the  run  can  be  restarted  simply  by  resubmitting  the 
original  command  file,  after  the  error  has  bean  corrected.  The  program  checks  for 
the  existence  of  the  statue  file  end  continues  execution  at  the  first  entry  which 
does  not  have  a date  and  time  entry,  whan  ell  paths  have  baen  processed,  the  program 
aalaulatee  the  field  strength  along  aach  path  using  the  parameters  specified  for  the 
transmitting  antenna  and  the  receiver.  These  data  are  written  to  a file  named  with 
the  extension  LWF  and  tha  status  file  ie  updated.  If  the  propagation  paths  were  eat 
up  automatically  by  ueer  specification  of  an  operating  area,  the  program  uses  the 
deta  in  the  LWF  fila  to  gansrate  a file  named  with  the  extension  CRD  and  the  status 
fila  ia  updated.  The  CRD  file  contains  values  of  the  signal  otrangth  and  its  stand- 
ard deviation  in  e grid  of  latitude  ve.  longitude  which  covers  the  operating  ares. 
This  last  file  may  now  be  used  in  a program  named  PL0T_CRD  to  obtain  geographical 
displays  of  the  signal  lavals,  signal  to  noisa  levels,  or,  together  with  other  CRD 
files,  signal  to  jammer  levels. 
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OUTPUT  DATA 


An  important  improvement  in  LWPC-1  is  in  the  handling  of  data  files.  Compared 
to  LWPC-0,  more  information  is  stored  in  these  files  and  the  formats  have  bean 
developed  to  allow  for  future  enhancements.  The  new  files  are  all  written  in  unfor- 
matted form  and  a epeoial  program  named  SCAN  hae  been  written  to  print  out  summarise 
of  the  parameters  stored  in  them.  Some  parameters  share  common  usage  throughout  the 
full  aet  of  output  data  files,  whenever  appropriate,  these  parameters  from  one  file 
are  passed  on  to  subsequent  types  of  files  in  order  to  provide  continuity  and  some 
audit  trail  information.  The  first  record  of  all  date  file  typae  contains  the  same 
information.  This  first  record  contains  an  8 oharacter  string  to  bs  used  to  rscord 
archive  information.  It  also  contains  three  stringe,  named  MDB-FILE»lrJ,  LWF-FILE-ID 
and  grd-FILE-id.  Each  of  these  strings  contains  the  date  the  file  was  written,  a 
randomly  generated,  3 oharacter  string  to  uniquely  identify  the  file  and  the  full 
file  name  including  the  directory  tree.  The  presence  of  these  strings  in  each  typo 
of  file  provides  information  regarding  the  history  of  the  data  whioh  was  usad  to 
produce  the  file.  For  inetanoe,  each  GRD  file  identifies  the  LWF  file  which  provided 
the  signal  strength  data  and  the  MDS  file  whiah  provided  the  mode  parameter  data 
used  to  calculate  the  LWF  data.  Generally,  the  program  identification  will  be  'LWPM- 
V10 1 . This  first  record  also  contains  the  list  of  paramsters  to  identify  the  propa- 
gation paths  for  the  data  aeta. 

A utility  program  named  SCAN  ia  used  to  produce  summaries  of  the  data  in  any 
of  the  output  files.  Examples  of  ths  output  from  this  program  using  the  sample 
problem  are  shown  below. 


SAMPLE  PROBLEM 

To  illustrate  the  capabilities  of  LWPC-1  a sample  case  is  presented,  The  run 
stream  is  shown  in  Figure  1.  An  operating  area  named  "Mediterranean"  is  used  to 
define  the  propagation  paths.  The  root  file  name  is  "SAMPLE"  so  that  the  following 
files  were  oreatadi  SAMPLE. eta,  SAMPLE. mds,  SAMPLE. lwf  and  SAMPLE jneditarranean.grd. 
The  tranemitter  parameters  were  retrieved  from  a specification -fi.le.  The  preview 
option  wee  exercised  to  produce  the  plot  shown  in  Figure  2.  A preview  option  was 
then  "turned  off"  by  indenting  the  aontrol  string  one  space. 

The  second  execution  of  the  program  produced  the  MDS,  LWF  and  GRD  files  named 
above.  Figure  3 shows  the  output  from  the  SCAN  program  for  the  first  path  in  the  MDS 
file.  The  output  identifies  ths  program  which  waa  used  to  generate  the  file  and  the 
date  that  it  was  run.  The  data  under  the  heading  which  begins  with  "mds"  is  a sum- 
mary of  the  parameters  for  each  seqment  along  the  path.  Figure  4 shows  the  output 
from  the  SCAN  program  for  the  LWF  file.  The  parameter  under  the  heading  "no"  identi- 
fies tho  number  of  field  components  in  the  mode  sums.  The  data  under  the  heading 
which  etarta  with  "bearng"  identify  the  pathu  by  their  bearing  angle  and  shows  the 
length  of  aacih  path.  The  signal  strength  as  a function  of  distance  from  the  trans- 
mitter along  the  northernmost  (bearing  - 24  degrees;  upper  panel)  and  southernmost 
(bearing  - 72  degrees;  lower  panel)  paths  are  plotted  in  Figure  5.  Along  the  bottom 
of  eaah  graph  is  a summary  of  the  important  path  segmentation  data,  namely,  the 
height  of  the  ionoophers  and  the  ground  conductivity.  The  beginning  of  each  segment 
is  indicated  by  a email  diamond  in  the  curve  representing  the  ground  conductivity. 
In  the  northern  path,  the  rapid  change  in  signal  strength  that  occurs  near  4 Mm  la 
due  to  the  e£  act  of  Greenland. 

Figure  o shows  the  SCAN  output  for  the  GRD  file.  This  output  identifies  the 
boundaries  of  tha  operating  area.  The  values  under  the  headings  "niat"  and  "nlon" 
are  the  number  of  latitudes  and  longitudes  used  to  define  the  grid.  Sample  contour 
maps  generated  by  PLOT  GRD  using  the  run  stream  in  Figure  7 era  presented  in  Figures 
e end  9.  The  first  of  £tiese  figure*  shows  contours  of  signal  strength  and  the  second 
shows  contours  of  signal  to  atmospheric  noine  ratio.  The  rapid  drop  in  signal 
strength  a»  paths  arose  Greenland  is  seen  in  the  closely  spaced  and  nearly  vertical 
contour  lines  in  the  eastern  pari  of  the  operating  area. 

DISTRIBUTION 

The  source  code  and  associated  data  sets  are  available  from  the  Defense  Tech- 
nical Information  Center. 
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$ run  Iwpa-dir; lwpm 

case-id  OMEGA  coverage  of  the  Mediterranean 

tx  SAMPLE 

tx-data  OMEGA-D  lwp: [dat] xmtr . list 

ionosphere  lwpm  day 

op-area  Mediterranean  30  10  45  -45 

lwf-vs-dist  ,11000,, 

print-swg  2 

print-lwf  2 

gcpath 
preview 

map-area  LANT  root  20  100  80  -50  7 5 

map-type  conductivity 

etart 

quit 


Figure  1;  Sample  run  stream  for  LWPM 
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Figure  3:  SCAN  Output  for  MDS  files 
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